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SUMMARY 


This  document  describes  the  package  of  special-purpose  computer  soft- 
ware delivered  and  installed  by  Goodyear  Aerospace  Corporation  (GAC)  in 
the  IJ.  Army  Engineering  Topographic  laboratories  (USA-ETL)  under 
Contract  DAAK70-76-C-02i»7. 

This  software  package  applies  an  area  correlation  function  to  two  in- 
put images: 

(1)  A live  image  taken  from  flight  tests  of  the  U.  S.  Army 
Pershing  II  Missile  Terminal  Guidance  System. 

(2)  A reference  image  for  use  with  the  above  guidance  system. 

The  software  package  has  been  operating  at  GAC  on  a HlS^/Xerox  Sigma 
9 computer  and  STARAN'’  parallel  processor.  The  subject  contract  provided 
for: 


(1)  Conversion  of  all  Sigma  9 assembly  language  code  and 
unacceptable  FORTRAN  code  to  FORTRAN  code  acceptable 
to  the  E.T.L.  Control  Data  Corporation  (CDC)  6ii00 
computer  system. 

(2)  The  conversion  of  all  operator  interactive  portions  of 

the  package  to  run  on  the  E.T.L.  DIAL  image  analysis  system. 

(3)  Conversion  of  appropriate  sections  of  code  to  that  applicabie 
to  the  E.T.L.  CDC  6^00/STARAN  interface. 

(A)  Preparation  of  engineering  documentation. 

(5)  Installation  and  demonstration  of  the  software  package 
at  E.T.L. 

(6)  Instructions  in  the  use  of  the  package. 

The  E.T.L.  Project  Monitor  for  this  contract  was  A.  T.  Blackburn, 
Advanced  Technology  Branch  of  the  Computer  Science  Laboratory.  The  GAC 
Project  Engineer  was  T.  E.  Gorsica,  Department  A75. 

NOTES: 

(1)  The  Program  Module  Descriptions,  Flow  Charts,  and  the  Program 
Listings  were  submitted  as  part  of  the  requirements  of  the 
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SECTION  I 
INTRODUCTION 


Correlation  guidance  systems  that  operate  on  the  principle  of  auto- 
matically comparing  a stored  reference  to  a sensed  image  of  a target  area 
.to  deter-r.ine  position  relative  to  the  target  offer  a number  of  advantages 
over  other  available  guidance  systems.  They  are  self  contained,  automatic, 
and  accurate  and  require  target  location  accuracy  only  relative  to  the 
immediate  target  area  without  requiring  a high  launch  poi nt-to-target 
accuracy.  A number  of  such  systems  are  being  developed  toward  operational 
deployment.  A common  requirement  for  all  of  these  systems  is  the  need  of 
a reference  for  each  target. 

The  ability  to  support  a guidance  system  operationally  with  suitable 
references  could  well  be  the  deciding  factor  on  the  effectiveness  of  these 
systems.  The  data  base  will,  in  turn,  be  the  key  to  the  reference  support. 
One  of  the  major  problems  in  the  guidance  reference  support  area  is  the 
establishment  of  requirements  for  a data  base.  The  data  base  materials 
which  have  in  the  past  been  used  would  be  expensive  to  prepare  and  use  for 
operational  support.  Little  consideration  was  given  to  operational  require- 
ments since  the  primary  objective  of  these  test  programs  was  to  demonstrate 
system  performance.  The  data  base  requirements  established  for  these  test 
programs  were  overly  stringent,  and  further  work  is  required  to  determine 
what  reduction  in  these  requirements  is  practical. 

The  program  provided  establishes  at  ETL  a capability  of  analyzing  the 
effect  of  introducing  reference  image  modification  in  terms  of  the  result- 
ant change  which  these  modifications  cause  in  the  live  versus  reference 
correlation  function.  This  capability  was  introduced  through  the  trans- 
ferral  of  an  existing  GAC*  software  system,  the  Precision  STARAN  Cor- 
relator. The  implementation  of  the  correlator  at  ETL  is  for  the  RADar 
Aimpoint  £uidance  (RADAG)®  system  since  there  is  an  urgent  requirement  for 
a definition  of  the  data  base  for  this  system.  However,  the  procedures 
can  be  readily  adapted  to  other  guidance  systems. 


^TM,  Goodyear  Aerospace  Corporation,  Akron,  Ohio  A1«3I5 
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SECTION  II 
METHOD  OF  SOLUTION 


The  basic  purpose  of  the  correlation  program  is  to  generate  cor- 
relation and  aperture  functions  of  digitized  imagery.  The  correlation 
function  is  given  by: 


’l.J 


N.N 


i.j-  -N,-N 


L.  .H.  •R..|  •«_! 
Ij  Ij  1+1, j+J 


(1) 


whe  re : 


M. 


i + l,j+J 
K 


is  a live- image  element 

is  a mask  eiement  (I  or  0) 

is  a reference -image  element 

is  the  number  of  elements  correlated 


The  aperture  function  required  to  obtain  normalized  correlation  is  given 
by: 

N.N 


c „ » 
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(2) 


where  C is  calculated  as  shown  below. 


A description  of  the  subscripts  used  in  Equation  (1)  and  (2)  is 
given  in  Figure  I.  The  vaiue  of  K is  obtained  using  the  relationship 


K - ir(R^  - rJ) 


(3) 


where  R2  is  the  outside  radius  of  the  mask, 


R|  is  the  inside  radius  of  the  mask. 


Correlation  and  aperture  calculations  can  be  conducted  along  any  row 
(generally  meaning  west  to  east)  or  column  (south  to  north).  The  value 
of  N can  be  selected  from  1 to  M.  Excursions  of  I and  J from  zero  are 
limited  by  the  size  of  the  live  scene  to  that  of  the  reference,  i.e., 
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I and  J cannot  exceed  plus  or  minus  (M-N) . This  limit  produces  a maximum 
of  2(M-N)+I  correlation  and  aperture  values  per  line.  The  mask  dimensions 
are  selectable  from  the  control  terminal. 

Computed  values  of  , and  Aj  j are  printed.  Additionally,  control 
parameters  and  i nput-image 'li ties  afe  printed. 

There  are  two  modes  of  operation  in  the  correlation  program,  viz., 
the  linear  mode  and  the  non-linear  (RADAG)  mode. 


a.  Linear  Mode  - The  linear  mode  uses  the  digitized  imagery  as  it  is 
received,  i.e.,  image  intensity  values  are  retained  throughout  the  process. 
Equations  (I)  and  (2)  are  used  in  this  mode.  The  value  of  C in  Equation 
(2)  is  the  maximum  value  of  the  live  image  within  the  correlation  aperture 
(defined  by  the  mask),  i.e., 


C 


Max 
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(M 


where  L and  M are  the  live  and  aperture  (masked)  data  sets,  respect i ve iy . 

b.  Non-Linear  (RADAG)  Mode  - The  non-linear  mode  of  operation  is 
programmed  to  simulate  an  adaptive  AGC  video  processor  and  a non-linear 
sample  weighting  process.  The  video  processor  measurements  are  simulated 
with  the  following  equations: 
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where: 

L ■ average  of  live  scene 
O ■ estimate  of  the  standard  deviation 


(6) 


The  effect  of  video  processor  AGC  and  the  non-linear  sample  weighting 
process  is  simulated  by  modifying  the  live-image  elements  with  the 
following  set  of  rules: 
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(7) 


IF  L. . > L + 2d 
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where  L'..  are  the  modified  live  image  elements. 


The  modified  live  image,  described  by  L',  is  subsequently  used  in 
Equation  (l)  to  compute  the  correlation  function.  Also,  in  the  non-linear 
mode,  the  value  of  C in  Equation  (2)  is  given  by: 

C - 2X0  . (8) 


This  factor  causes  the  quotient  ^/A,  designated  in  the  next  section, 
to  be  proportional  to  the  percentage  of  data  correlating  at  any  point. 

c.  Output  Data  - A plot  program  uses  the  j and  Aj  j values  from 
the  output  of  the  correlation  program.  For  each  line  processed, 
the  normalized  correlation  function  i{/|  j is  computed  using 


(9) 


Additionally,  the  derivatives  of  lli,  namely,  ip'  and  iK'  are  computed  along 
the  line  chosen.  The  averages  and  standard  deviations  of  (J,  A,  and  i|), 
are  computed  and  printed  out  along  with  the  live  and  reference  image  titles. 
The  maximum  and  minimum  values  of  A,  i|i'  and  ip”  are  printed  out  as 
well  as  the  positions  of  the  maximum  and  minimum  values  of  ip,  ip',  and  ip". 
Additionally,  a three-point  parabolic  curve  fit  program  has  been  provided 
to  predict  the  locations  of  the  maximum  and  minimum  values  of  ip  with  a 
positional  resolution  better  than  one  integer  pixel. 


SECTION  I 1 1 
RESOURCES  REQUIRED 


The  GYPSC  program  opeates  under  the  ETL/OIAL  interactive  environment- 
The  FORTRAN  program  utilizes  the  COC  6400,  GAC  STARAN,  and  PDP/COMTAL/ 
TEKSCOPE  computer  systems  in  concert  with  their  established  ETL  facility 
hardware  and  software  interface  subsystems. 


The 

user  application  program  GYPSC  includes  the  MAIN  program  of  the 

same  name  and  i ts 

21  associated  subroutines  as  follows: 

1. 

GYPSC 

12. 

CORRELV 

2. 

YESNO 

13. 

FETCH 

3. 

GETINP 

14. 

GETREF 

4. 

DECODE 

15. 

SWAP 

5. 

BIGTOP 

16. 

TURN 

6. 

OMENU 

17. 

SKIPR 

7. 

TEKOUT 

18. 

LIVBLD 

B. 

GETSCN 

19. 

IRADAG 

9. 

STAT 

20. 

RADAGF 

10. 

PLOTR 

21. 

APERTUR 

11. 

CORREL 

22. 

APSET 

Fi rst-level 

FORTRAN/SYSTEM/ library  calls  include: 

I. 

CLOSMS 

7. 

EOR 

2, 

READMS 

8. 

OPENMS 

3. 

MOVLEV 

9. 

REWIND 

4. 

READEC 

10. 

WRITMS 

5. 

WRITEC 

6. 

ENCODE 

Fi rst-level 

DIAL-FORTRAN  interface  subroutine  calls  include: 

1. 

FINIS 

8. 

DREAD 

2. 

LBLRD 

9. 

FIND 

3. 

LOCATE 

10. 

GSET 

4. 

PMINIT 

11. 

IMGDSK 

5. 

TEKMSG 

12. 

LBLEXT 

6. 

TITLE 

13. 

SETCORE 

7. 

UNPKI 

14. 

TEKRD 

The  Q9PFMGR  routine  is  also  called  at  first-level  to  link  the  perma- 
nent LUN  (Logical  Unit  Number)  F:94  FORTRAN  file  which  stores  the  menu 
fi le.  MENU. 


First-level  STARAN- INTERFACE  software  (refer  to  GER-16352)  calls  in- 
clude: 


1. 

SCNTRL 

2. 

RCNTRL 

3. 

STATUS 

1*. 

LOADPRG 

5. 

DETACH 

As  configured,  the  GYPSC  program  requires  a CM  (central  memory) 
execution  field  length  (FL)  of  106,000  (octal).  The  ECS  (extended  core 
storage)  requirement  is  36,861*  (decimal).  Two  (2)  COMTAL  CRT  display 
units  are  used. 


SECTION  IV 
OPERATION 


6YPSC  is  run  as  an  operator  interactive  program  from  the  DIAL  system 
terminal.  Before  loging  on  the  terminal,  the  operator  must  insure  that 
the  STARAN  and  the  CDC  61*00  are  ready  for  communication  with  one  another. 
This  is  accomplished  by: 

(1)  Turning  "ON"  CDC  peripheral  device  #31 
(STARAN)  via  the  CDC  E display; 

(2)  Invoking  the  CDC  communications  program  in  the 
ST  .ARAN  by  typing  "BA  CDC  10  (CR)"  after  the 

$ (dollar  sign)  on  the  STARAN  decwriter. 

GYPSC  is  resident  on  the  CDC  disk  memory  as  an  absolute  program.  It 
runs  in  a field  length  of  106,000  (octal)  words,  less  than  half  of  the 
available  CDC  core  memory.  GYPSC  is  accessed  by  loging  on  to  DIAL  in  the 
normal  way,  selecting  major  projects  sub-element  6 (PROGRAM  DEVELOPMENT), 
and  typing  in  the  program  name,  "GYPSC",  as  illustrated  in  Figure  2. 

Should  GYPSC  no  longer  be  resident  on  the  CDC  disk,  it  is  only  necessary 
to  submit  the  CDC  program  card  deck  for  compilation.  All  control  cards 
are  included  in  the  deck. 

GYPSC  asks  for  the  DIAL  file  image  names  of  two  images  to  be  cor- 
related, as  shown  in  Figure  3.  Note  that  figures  referenced  in  this 
section  are  reproductions  of  actual  GYPSC  runs  taken  from  the  DIAL  system 
TEKSCOPE.  Operator  entries  have  been  underlined  on  these  reproductions  to 
distinguish  them  from  program  generated  text.  Should  the  operator  enter 
a non-existent  image  name,  he  will  be  notified  and  GYPSC  will  again  ask 
for  a valid  name.  (Should  a CR  be  entered  instead  of  a valid  name,  the 
DIAL  system  software  will  present  a catalogue  of  valid  DIAL  images.  The 
system  software  will  return  to  GYPSC  if  a CR  is  entered  into  the  catalogue 
program.)  When  the  two  valid  entries  are  made,  GYPSC  responds  by  placing 
the  LIVE  image  on  COMTAL  A,  the  REFERENCE  image  on  COMTAL  B,  displaying 
the  image  header  information  (if  any)  and  asking  for  further  directives. 

Note  that  the  images  used,  GYLIVE  and  GYREF,  did  not  have  any  header 
information  (Figure  21  shows  two  images  that  contain  header  information). 
GYLIVE  and  GYREF  are  identical  images  except  that  the  pixel  intensities 
are  reversed  (one  image  is  the  negative  of  the  other). 

The  "MENU"  directive  calls  up  a listing  of  the  parameters  and 
directives  available.  The  "SUMM"  directive  displays  the  current  settings 
of  all  of  the  GYPSC  input  PARAMETERS.  Figure  4 shows  the  MENU  and  a SUMMARY 
of  the  default  input  parameters  which  are  the  values  in  effect  when  GYPSC 
is  first  accessed,  and  until  changed  by  subsequent  operator  directives. 

The  run  number  is  incremented  each  time  the  "START"  directive  is 
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entered  to  begin  the  designated  correlation(s) . 

A "CROSS"  correlation  performs  an  area  correlation  of  a subset  of  the 
LIVE  image  against  the  '^hole  of  the  REFERENCE  image.  An  "AUTO"  correlation 
performs  an  area  correlation  of  a subset  of  the  REFERENCE  image  against 
the  whole  of  the  REFERENCE  image. 

The  axis  may  be  HORIZONTAL,  VERTICAL,  or  BOTH.  A HORIZONTAL  cor- 
relation selects  a doughnut  shaped  patch  (designated  the  aperture)  of  the 
LIVE  image  with  its  center  at  L ,L  and  performs  an  area  correlation  of 
this  patch  against  all  of  the  possible  areas  of  the  REFERENCE  image  that 
Tall  on  a horizontal  line  passing  through  R ,R  on  the  REFERENCE  image. 

See  Figures  5 and  6.  Stated  another  way,  tfie  J(atch  of  the  LIVE  image  is 
slid  horizontally  a pixel  at-a-time  across  the  REFERENCE  image  from  west 
to  east.  At  each  possible  place  the  patch  from  the  LIVE  image  is  area 
correlated  against  the  underlying  area  of  the  REFERENCE  image.  This 
generates  a curve  of  correlation  amplitude  versus  the  horizontal  position 
of  the  patch  of  the  LIVE  image.  A VERTICAL  correlation  similarly  slides 
the  patch  of  the  LIVE  image  vertically  from  the  south  to  the  north  passing 
through  point  R ,R  on  the  REFERENCE  image.  When  "BOTH"  is  entered, 
both  a HORIZONTAL  Xnd  VERTICAL  correlation  are  successively  performed. 

In  the  RADAG  mode,  the  Intensities  of  the  LIVE  image  are  re-mapped 
by  an  algorithm  that  simulates  the  adaptive  A6C  video  processor  and  a non- 
linear sample  weighting  process  in  the  hardware  RADAG  correlator.  In  the 
LINEAR  mode,  no  change  is  made  to  the  intensities  of  the  LIVE  image. 

The  number  of  lines  of  video  sent  to  the  STARAN  processor  is  deter- 
mined by  "LINES".  LINES  should  be  equal  to  or  greater  than  R^,  lest  the 
edges  of  the  round  doughnut  shaped  aperture  become  flattened.  However, 
if  LINES  is  set  to  an  overly  large  value,  the  program  execution  time  will 
increase  unnecessarily. 

When  the  LIVE  image  comes  from  the  guidance  system's  radar,  (the 
usual  operating  mode),  Rj  can  be  set  to  the  inner  radius  of  the  scene 
content  (the  minimum-ground-range  features  picked  up  by  the  radar).  R2 
can  be  set  to  the  maximum  ground  range.  All  magnetic  tapes  of  live 
radar  imagery  generated  on  the  Goodyear  Aerospace  radar  re-sampler  system 
have  a property  that  can  be  used  to  automatically  set  the  aperture  in 
GYPSC.  These  tapes  have  each  pixel  coded  as  an  8-bit  byte.  When  the 
resampler  determines  that  a pixel  is  not  between  the  minimum  and  maximum 
radar  ground  range,  the  pixel  value  is  set  to  zero.  When  the  value  of 
R-'R.  is  equal  to  150,  the  GYPSC  program  will  treat  a live  pixel  as  being 
within  the  aperture  if  it  has  a non-zero  intensity  value.  This  lets  the 
scene  content  of  the  resampled  live  image  accurately  set  the  aperture  size. 
This  also  implies  that  there  should  be  no  zero  valued  pixels  within  the 
aperture  of  the  live  images.  As  a special  case,  if  the  difference  between 
R-  and  R^  is  equal  to  150,  GYPSC  will  use  a square  shaped  aperture  of  301 
pixels  on  a side.  This  may  be  found  useful  when  correlating  full  area 
images  not  pertinent  to  the  RADAG  guidance  system.  Otherwise  Rj  and  R^ 
prescribe  the  circular  aperture. 

The  lambda  parameter  is  the  integer  number  of  standard  deviations 
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accepted  in  the  RAOAG  mode.  See  Equation  (7)  in  the  previous  section. 

Figure  7 shows  entries  that  change  R2 > LINES  and  AXIS  followed  by 
another  summary.  Param;;ters  may  be  entered  in  any  order  and  may  be  re- 
entered any  number  of  times.  They  will  remain  constant  from  run-to-run 
until  changed  by  the  operator.  The  START  command  initiates  the  actual  cor- 
relation, the  results  of  which  are  shown  in  Figure  8.  Note  that  each  time 
GYPSC  erases  the  TEKSCOPE  screen,  the  ledger  is  reprinted  at  the  top  so 
that  all  of  the  pertinent  information  appears  along  with  the  correlation 
resul ts. 

Figure  8 also  shows  operator  entries  for  width,  offset,  and  curve 
number.  The  curve  of  correlation  amplitude  returned  from  the  STARAN  pro- 
cessor, shown  in  Figure  9A,  contains  information  over  a greater  excursion 
than  is  typically  of  interest.  The  curve  displayed  on  the  TEKSCOPE  is  a 
subset  of  the  total  correlation  curve  as  shown  in  Figure  98.  The  subset 
is  selected  by  the  width  and  offset  entries.  The  TEKSCOPE  screen  size 
limits  the  "WIDTH"  entry  to  a maximum  value  of  5^.  This  subset,  or  window, 
also  describes  the  limits  of  the  correlation  curve  over  which  numerical 
statistics  of  the  correlation  process  are  taken.  The  "CURVE  NUMBER"  entry 
describes  which  of  five  possible  types  of  curves  (as  defined  in  the  menu 
shown  in  Figure  are  plotted  on  the  TEKSCOPE.  The  statistics  include 
the  maximum,  minimum,  average,  and  standard  deviation  of  the  unnormalized 
correlation  function,  the  aperture  function,  and  the  normalized  correla- 
tion function  (PSI).  The  statistics  also  include  the  minimum  and  maximum 
values  of  the  first  derivative  of  PSI  (DPS I)  and  the  second  derivative  of 
PSI  (DDPSI),  as  well  as  the  positions  of  the  minimum  and  maximum  values  of 
PSI,  DPSI,  and  DDPSI.  The  statistics  are  dependent  on  the  window  selected 
by  the  WIDTH  and  OFFSET  but  do  not  depend  on  the  curve  number  entered. 

Note  that  there  is  a local  minimum  near  the  center  in  the  plot  of  the 
un-normal i zed  correlation  function  (see  Figure  8).  The  center  of  the  plot 
is  always  denoted  by  a plus  (+)  sign.  Figures  9A  and  98  tie  the  center  of 
the  plot  back  to  the  origin  of  the  reference  image  (Figure  5)  positionally. 
As  the  two  example  images  (GYLIVE  and  GYREF)  are  of  opposite  polarity,  the 
local  minimum  indicates  the  lateral  position  at  which  the  two  images  are 
most  dissimilar.  (This  is  also  known  as  the  quotient  match  point.)  The 
minimum  is  only  local  because  the  un-normal i zed  correlation  function  is 
scene  dependent.  If  the  operator  responds  to  the  question  "SHALL  I 
CONTINUE"?  with  a carriage  return,  as  was  done  in  Figure  8,  GYPSC  will 
give  the  operator  the  opportunity  to  re-enter  the  statistics  without  re- 
running the  whole  correlation  program.  Figure  10  shows  the  aperture 
function  (curve  number  2)  for  the  same  run.  Looping  back  to  the  statistics 
again  and  entering  curve  number  3 yields  a plot  of  the  normalized  cor- 
relation function  PSI  (Figure  11).  Note  that  the  minimum  is  much  more 
pronounced  because  the  gross  scene  dependencies  have  been  removed  in  the 
normalization  process.  The  strength  of  the  match  is  indicated  by  the 
I difference  between  the  minimum  and  maximum  values  of  PSI.  The  position 

► of  the  match  is  indicated  by  the  position  of  the  minimum  value  of  PSI 

which  is  reported  as  plus  I pixel  on  an  integer  basis  and  as  plus  0.28 
pixel  by  a three  point  parabolic  curve  fit  interpolator. 

Figure  12  shows  curve  number  k,  the  first  derivative  of  PSI,  and 
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Figure  13  shows  curve  number  5>  the  second  derivative  of  PSI.  Figure  lA 
shows  the  plot  of  PSI  obtained  by  entering  an  offset  other  than  zero  while 
Figure  15  shows  PSI  with  a larger  width.  Note  that  the  statistics  reported 
in  these  last  two  figures  have  changed  because  the  size  and  position  of 
the  window  have  changed.  The  window  may  be  enlarged  to  show  the  activity 
on  the  skirts  of  the  curve,  or  restricted  to  confine  the  statistics  to  a 
local  area  of  the  curve.  At  the  bottom  of  Figure  15  the  operator  entered 
"Y"  to  continue  instead  of  looping  back  to  re-enter  the  statistics. 

Figure  16  shows  the  vertical  normalized  correlation  curve  of  the  same 
run.  The  operator  may,  at  the  bottom,  enter  a carriage  return  to  loop 
back  and  re-enter  the  statistics.  However  the  entry  of  a "Y"  or  a "N" 
will  take  him  back  to  the  top  of  the  GYPSC  program  because  Run  1 is 
complete. 

Figure  17  shows  the  beginning  of  the  GYPSC  program.  Since  GYPSC  knows 
that  it  has  already  been  run  once,  it  gives  the  operator  the  option  of 
terminating  completely,  using  the  images  already  inserted,  or  respecifying 
the  images.  In  this  example  the  operator  chose  the  latter.  By  entering 
"AI"  for  COMTAL  A or  "Bl"  for  COMTAL  B,  either  image  already  existent  on 
either  COMTAL  may  be  used  for  the  live  image,  or  for  the  reference  image. 
This  example  then  shows  that  GYPSC  was  directed  to  run  an  auto-correlation 
in  the  horizontal  direction  only  and  to  give  a summary  of  the  parameters 
before  actually  starting  the  run. 

Because  an  auto-correlation  is  a correlation  of  a portion  of  the 
reference  image  against  the  whole  reference  image,  the  resultant  normal- 
ized correlation  function  as  shown  in  Figure  18  gives  a peak  at  the 
position  of  best  correlation,  plus  0.12  pixels  in  this  example. 

Figure  19  directs  a RADAG  cross-correlation.  The  results  shown  in 
Figure  20  are  interpreted  the  same  way  as  the  linear  mode  except  that 
additional  statistics  about  the  live  scene  are  given.  AVEL  is  the  average 
value  of  the  original  live  scene  pixels  that  are  within  the  aperture. 

SIGMA  is  an  estimate  of  the  standard  deviation  of  these  live  scene  pixels. 
RMAXL  equals  AVEL  + ( LAMBDA- S I GMA) . RMINL  equals  AVEL  - ( LAMBDA- S 1 GMA) . 
LMINis  the  minimum  value  of  these  live  scene  pixels  and  LMAX  is  the  maximum 
value  of  these  live  scene  pixels. 

Figure  21  shows  how  the  proper  parameters  are  taken  from  the  header 
information  encoded  on  tapes  of  live  radar  imagery  generated  on  the 
Goodyear  Aerospace  radar  re-sampler  system.  The  normal  use  of  the  GYPSC 
program  is  to  correlate  a re-sampled  live  radar  image  against  a reference 
image.  Because  an  appropriate  pair  of  resampled  live  and  reference  images 
were  not  available  in  DIAL  image  format  when  these  documentation  runs  were 
made,  two  re-sampled  live  images  were  used,  (PPIGYl)  and  (PPIGY2). 

These  two  images  were  taken  on  the  same  run  of  the  same  flight  but  differ 
slightly  in  altitude.  See  Table  I for  definitions  of  the  legends  used  in 
live  radar  image  headers. 


r 


-21- 


GVPSC  •</N/77  it.te.K 


SMALL  1 CONTHAJC^Oyt  M tHt|  LOOM  OM  CR  ONLVI 


GYPSC  •4/29/77  ll.t3.42. 

STARAN/CDC64de  AREA  CORRELATION 

.IJIC  IMCC CVLIVC 


I 

.j 


SK 

I 


># 

I 


>m 

I 


u. 

25 

XO 

«« 

o 

• 

• 

• 

VI 

ifS 

i 

>ft 

KO 

U» 

if 

5 

oc 

K 

! 


U. 

5 


VI 

K 


I 


IA»«» 


8 

oyi<Ap> 
m •tfi 

vif*i«*  • 


ftC 

I? 


VI 

'lAIOl 


C ••noiA 
S#  •Vi/kVI 
C •r* oi^f^ 
•^CBAI^OIOI 

« axAnoi  •« 
C V — • . 4 


8S 

t I 


Aim  mA-  tt 

mm 

II  *11 

8 


I 


iifKf 

Ul 

1 

iz 

4A 

C •ao««« 

«-»4^AimAin 


mm^ 

f-  A#^ 

• t I 


01 

I 


lAVIVk 


2S 


I 


. iA  VI VIIA 

“-sa 

O VI  o 
OVI 
0.0^ 
a.Q 


iH^tL  I C0NTlNUC7(tV>  M MOi  tW  OM  « ONLY) 


GYPSC  •4ya»/77  I1.14.M. 

STARAN/CDC64ee  AREA  CORRELATION 


Figure  15 


Figure  17 


GVPSC  94/29/77  It.Il.M 


« 

Ui  w 

|V  9 

“5  I 

O • 

A 


I 


Ain<«r^ 


ui^ 

Iss". 


<VAi  #r»*« 


« lAl  ▼ 

8 i |882«5  ? 

zm  ► wn  m 


• * • f 


8S  S72 

SS82^!S 

C • 

«V 

S^OArw^J  «« 

Kiwitf  • • • 

I I 


Su#***"**  ■••••• 

OavA«/>V) 

r?‘8&  “•&& 

«B  O ^ A 

pl#» 


SIMLl  I C0HTlhUC7(lVt  M *m|  LOO^  OM  Ca  ONLY  I 


tSTMT 


TABLE  I.  RADAR  IMAGE  HEADER  LEGENDS 


IMAGE  HEADER 

LEGEND 

FLT 

FI  ight  number 

RUN 

Run  number 

CPI 

Checkpoint  initiation 

VT 

Time  and  date  when  the  video  tape 
was  digitized 

PR 

Time  and  date  when  the  digitized 
image  was  converted  from  polar  to 
rectangular  coordinates 

FLT  TIME 

IRIG  standard  time  off  of  the  video 
tape 

CPI  1 

Initial  CPI 

CPIF 

Final  CPI 

Al 

Initial  altitude  (meters  A.G.L.) 

AF 

Final  altitude  (meters  A.G.L.) 

Nl 

Initial  position  of  nadir  point 
north  of  the  target  (meters) 

NF 

Final  position  of  nadir  point  north 
of  the  target  (meters) 

El 

Initial  position  of  nadir  point 
east  of  the  target  (meters) 

EF 

Final  position  of  nadir  point  east 
of  the  target  (meters) 

ANGI 

Initial  azimuth  angle  (degrees  w.r.t. 
true  north)  of  the  antenna  n*ien 

CPI  started 

B 

Band  number 

L is  obtained  by  divding  east  initial  (El)  by  the  scale  factor  (in 
meters/^pixel) . L is  obtained  by  dividing  north  initial  (Nl)  by  the 
scale  factor.  Thi(  scale  factor  will  be  a constant  for  each  band  and 
should  be  requested  when  re-sampled  live  images  are  obtained. 

R and  R can  initially  be  made  the  sar^  as  and  L , respectively. 
If  an  fni t ial^vert i cal  correlation  shows  an  appreciable  positional  error, 

R can  be  adjusted  correspondingly  for  a subsequent  horizontal  correla- 
tYon  to  obtain  a stronger  horizontal  correlation  peak.  Conversely,  if  an 
initial  horizontal  correlation  shows  an  appreciable  position  error,  R^  can 
be  adjusted  correspondingly  for  a stronger  vertical  correlation  peak. 
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The  PSI  curves  generated  by  GYPSC  for  horizontal  and  vertical  correlations 
can  be  thought  of  as  cuts  through  the  correlation  surface. 

Note  that  due  to  the  nature  of  a digital  correlation  program  such  as 
GYPSC,  L and  L can  only  be  integer  numbers*  Additional  positional 
accuracy^may  be^obtained  by  modifying  the  correlation  position  reported 
by  GYPSC  by  the  remainders  obtained  from  dividing  El  and  Nl  by  the  scale 
factor. 

Figure  22  reports  that  the  portion  of  live  Image  PPIGYl  centered  at 
X-  - 1 correlates  best  at  position  X - -13.91.  ((POS  PSI  MAX) -OFFSET) , 

in  the  reference  image  PPIGY2. 

Figure  23  shows  one  method  of  terminal Ing  the  GYPSC  program.  The 
other  method  is  to  enter  "DONE"  when  GYPSC  is  requesting  commands  and 
di recti ves. 
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